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ABSTRACT 

The association of 7-ray bursts (GRBs) and core-collapse supernovae (SNe) 
of Type lb and Ic was motivated by the detection of SN 1998bw in the error box 
of GRB 980425 and the now-secure identification of a SN 1998bw-like event in the 
cosmological GRB 030329. The bright radio emission from SN 1998bw indicated 
that it possessed some of the unique attributes expected of GRBs, namely a 
large reservoir of energy in (mildly) relativistic ejecta and variable energy input. 
The two popular scenarios for the origin of SN 1998bw are a typical cosmological 
burst observed off-axis or a member of a new distinct class of supernova explosions 
(gSNe). In the former, about 0.5% of local Type Ib/c SNe are expected to be 
similar to SN 1998bw; for the latter no such constraint exists. Motivated thus, 
we began a systematic program of radio observations of most reported Type Ib/c 
SNe accessible to the Very Large Array. Of the 33 SNe observed from late 1999 to 
the end of 2002 at most one is as bright as SN 1998bw. From this we conclude that 
the incidence of such events is < 3%. Furthermore, analysis of the radio emission 
indicates that none of the observed SNe exhibit clear engine signatures. Finally, 
a comparison of the SN radio emission to that of GRB afterglows indicates that 
none of the SNe could have resulted from a typical GRB, independent of the initial 
jet orientation. Thus, while the nature of SN 1998bw remains an open question, 
there appears to be a clear dichotomy between the majority of hydrodynamic 
and engine-driven explosions. 

Subject headings: radio continuum:supernovae — supernovae:type Ib/c — gamma 
rays:bursts 
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1. Introduction 

The death of massive stars and the processes that lead to the formation of the compact 
remnants is a forefront area in stellar astrophysics. Recent advances in modeling suggest that 
great diversity can be expected. Indeed, observationally we have already witnessed a large 
diversity in the neutron star remnants: radio pulsars, AXPs, SGRs, and the central source 
in Gas A. We know relatively little about the formation of black holes, static or rotating. 

The compact object form following the collapse of the progenitor core. The energy of the 
resulting explosion could be supplemented or even dominated by the energy released from the 
compact object (e.g. a rapidly rotating magnetar or an accreting black hole). Such "engines" 
can give rise to asymmetrical explosions (MacFadyen & Woosley 1999). Regardless of the 
source of energy, a fraction of the total energy. Ex, is coupled to the debris or ejecta (mass 
Mgj) and it is these two gross parameters which determine the appearance and evolution of 
the resulting explosion. Equivalently one may consider Ek and the mean initial speed of 
ejecta, Vq, or the Lorentz factor, Fq = [1 — Po]~^^'^, where Po — Vq/c. 

Supernovae (SNe) and 7-ray bursts (GRBs), are distinguished by their ejecta velocities. 
In the former vq ~ 10^ km s~^ as inferred from optical absorption features (e.g. Filippenko 
1997), while for the latter Fq > 100, inferred from the non-thermal prompt emission (Good- 
man 1986; Paczynski 1986), respectively. The large difference in initial velocity arises from 
significantly different ejecta masses: Mej ~ few in SNe compared to ~ 10~^ in 
GRBs. 

In the conventional interpretation, Mpj for SNe is large because Ek is primarily derived 
from the (essentially) symmetrical collapse of the core and the energy thus couples to all 
the mass left after the formation of the compact object. Mysteriously, Ek clusters around 1 
FOE (FOE is 10^^ erg) in most SNe, a mere 1% of the energy released in the gravitational 
collapse of the core. 

Whereas the initial ejecta speed is solely determined by Ek and Mej, a fraction of the 
ejecta is accelerated to higher velocities as the blast wave races down the density gradient 
of the stellar enveloped (e.g. Matzner & McKee 1999). For the wind- or binary- stripped 
(e.g. Uomoto 1986; Branch, Nomoto & Fihppenko 1991; Woosley, Langer & Weaver 1993; 
Nomoto et al. 1994) hehum and carbon progenitors of Type lb and Ic SNe, both factors (a 
smaller core mass and a steep density gradient) conspire to produce ejecta at velocities as 
high as TjS ~ 1. However, only < lO^^E'x is carried by these ejecta (Golgate 1968; Woosley & 
Weaver 1986; Matzner & McKee 1999). In contrast, high velocity ejecta is neither expected 
nor observed in Type II SNe with their massive stellar envelopes. 

GRB models, on the other hand, appeal to a stellar mass black hole remnant, which 
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accretes matter on many dynamical timescales and powers relativistic jets (the so-called 
coUapsar model; Woosley 1993; MacFadyen & Woosley 1999); highly magnetized neutron 
stars have also been proposed (e.g. Ruderman, Tao & Kluzniak 2000). Observationally, this 
model is supported by the association of some GRBs with SN explosions (e.g. Stanek et al. 
2003). In addition, the complex temporal profiles and long duration of GRBs are interpreted 
in terms of an engine that is relatively long lived (i.e. not a singular explosion). The high 
Lorentz factors, a high degree of coUimation with opening angles of a few degrees (Frail et 
al. 2001), and episodes of energy addition presumably from shells of ejecta with varying 
Lorentz factors, further distinguish GRBs from Type Ib/c SNe. 

We now recognize that engine-driven events - GRBs and the recently discovered X-ray 
Flashes (Heise, in 't Zand & Kulkarni 2003) - in fact have a wide dispersion in their ultra- 
relativistic output as manifested by their beaming-corrected 7-ray energies (Bloom, Frail 
& Kulkarni 2003) and X-ray luminosities (Berger, Kulkarni & Frail 2003). However, these 
cosmological explosions appear to have a nearly constant total explosive yield when taking 
into account the energy in mildly relativistic ejecta (Berger et al. 2003). 

The unusual and nearby {d ~ 40 Mpc) SN 1998bw shares some of the unique attributes 
expected of GRBs. This Type Ic SN was found to be coincident in time and position with 
GRB 980425 (Galama et al. 1998), an event with a single smooth profile. The inferred 
isotropic energy in 7-rays of GRB 980425 was only 8 x 10^'' erg (Plan et al. 2000), three to 
six orders of magnitude fainter than typical GRBs. More importantly, SN 1998bw exhibited 
unusually bright radio emission indicating about 10^° erg of mildly relativistic ejecta (Li & 
Ghevalier 1999). Equally significant, the radio emission indicated a clear episode of energy 
addition (Li & Chevaher 1999). None of these features - 7-rays, significant energy with 
r/3 > 2, and episodes of energy addition - have been seen in any other nearby SN. Thus, 
the empirical data strongly favor an engine origin for SN 1998bw. In the optical, SN 1998bw 
also appears to be extreme: velocity widths approaching 60,000 km s~^ were seen at early 
time (Iwamoto et al. 1998) and the inferred explosion energy may be above normal values, 
with estimates ranging from 2 to 50 FOE (Hoflich, Wheeler & Wang 1999; Nakamura et al. 
2001). 

The inference of an engine in SN 1998bw raises two scenarios for its origin and relation 
to GRBs. GRB 980425 may have been a typical burst but viewed well away from the jet axis 
(hereafter, the off-axis model), thereby resulting in apparently weak 7-ray emission despite 
the great proximity. Alternatively, SN 1998bw represents a different class of SNe. If so, 
coUapsars can produce very diverse explosions. 

A powerful discriminant between these two scenarios is the expected rate of SN 1998bw- 
like events. In the off-axis model, the fraction of Type Ib/c SNe that are powered by a 
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central engine is linked to the mean beaming factor of GRBs (e.g. Frail et al. 2001; Totani 
& Panaitescu 2002). Recently, Frail et al. (2001) presented the distribution of jet opening 
angles for a sample of 15 GRBs, and found that the mean beaming factor is (/^~^) ~ 500; 
here fb = [1 — cos{6j)] is the beaming fraction, and 6j is the collimation angle. With an 
estimated local GRB rate of ~ 0.5 Gpc^'^ yr~^ (Schmidt 2001) compared to a Type Ib/c 
SN rate of ~ 4.8 x 10^ Gpc'^ yr^^ (Marzke et al. 1998; Cappellaro, Evans & Turatto 1999; 
Folkes et al. 1999), we expect that ~ 0.5% of IVpe Ib/c SNe will be similar^ to SN 1998bw. 

On the other hand, if SN 1998bw is not an off-axis cosmological burst, then the rate 
of similar events has to be assessed independent of the GRB rate. An upper limit can be 
obtained by assuming that all Type Ib/c SNe are engine driven highly asymmetric explosions 
with SN1998bw having the most favorable geometry. In this context, Norris (2002) has 
argued that of the 1429 long-duration BATSE bursts, about 90 events possess similar high- 
energy attributes as that of GRB 980425. This sub-sample may be concentrated along the 
super-galactic plane. If this sub-sample is accepted as distinct from the cosmological bursts 
then ~ 25% of Type Ib/c SNe within 100 Mpc are expected to be events like SN 1998bw. 

Here, we report a comprehensive program of radio observations of nearby Type Ib/c 
SNe. We began this program in 1999 (motivated by SN 1998bw) and observed most reported 
Type Ib/c SNe with the Very Large Array. Our basic hypothesis is that (mildly) relativistic 
ejecta are best probed by radio observations, as was demonstrated in the case of SN 1998bw. 
Furthermore, radio observations of Type Ib/c SNe allow us to directly compare these objects 
to the radio afterglows of cosmological GRBs. Thus, we can empirically (direct comparison 
of radio luminosity distributions) and quantitatively (calorimetry via radio observations) 
investigate the hnk, or lack thereof between Type Ib/c SNe and cosmological GRBs. As 
alluded to above, we did not investigate Type II SNe since the extended envelopes and dense 
circumstellar media of their progenitors are reasonably expected to mask the activity of a 
putative engine and thus suppress the presence of mildly relativistic ejecta to which we are 
most sensitive. 

The organization of the paper is as follows. In §2 we present the details of the obser- 
vations. The results are summarized in §3, where we investigate the broad radio properties 
(§3.1), expansion velocities (§3.2), and energies in high velocity ejecta (§3.3). We further 
provide a comparison to the radio afterglows of GRBs in §4 and draw conclusions in §5. 



■^We note that this fraction may be somewhat higher in jet models in which the energy and/or Lorentz 
factor decrease away from the jet axis. The exact fraction depends on the details of the energy and Lorentz 
factor distribution (Rossi, Lazzati, & Rees 2002; Zhang & Meszaros 2002). 
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2. Observations 

Tabic 1 summarizes the Very Large Array (VLA^) observations of Type Ib/c SNe start- 
ing in late 1999 and up to the end of 2002. We observed a total of 33 SNe out of 51 
identified spectroscopically during the same period. The observed targets were determined 
solely by the availability of observing time and optical selection criteria; we did not employ 
any additional selection criteria. 

In all observations we used the standard continuum mode with 2 x 50 MHz bands, 
centered on 1.43, 4.86, or 8.46 GHz. We used the sources 3C48 (J0137+331), 3C 147 
(J0542-I-498), and 3C286 (J1331-I-305) for flux cahbration, and calibrator sources within 
~ 5° of the SNe to monitor the phase. The data were reduced and analyzed using the 
Astronomical Image Processing System (Fomalont 1981). 

3. Population Statistics 

In this section we investigate the ejecta properties and diversity of the sample. Results 
for individual SNe are given in the Appendix. In Figure 1 we plot the radio luminosities 
and upper limits for Type Ib/c SNe observed in this survey and in the past (SN 1983N: 
Sramek, Panagia & Weiler 1984; SN 1984L: Panagia, Sramek & Weiler 1986; Weiler et al. 

1986; SN1990B: van Dyk et al. 1993; SN 1998bw: Kulkarni et al. 1998; SN2002ap: Berger, 
Kulkarni & Chevalier 2002). The typical delay between the SN explosion and time of our 
observations is about 20 days, with four SNe observed with a delay of over 100 days. In 
addition, three of the SNe are embedded in host galaxies with strong radio emission. For 
these cases, we adopt upper limits that correspond to the brightness of the galaxy (at least 
ten times the root-mean-square noise of the individual image) . Four of the thirty three SNe 
have been detected. Thus the detection rate of our experiment with a typical flux density 
limit of 0.15 mJy (3(t) is about 12%. 

3.1. Radio Properties of Type Ib/c SNe 

Figure 1 provides a succinct summary of the radio light curves of the Type Ib/c SNe. 
Two strong conclusions can be immediately drawn from this Figure. First, SNe as bright 



^The VLA is operated by the National Radio Astronomy Observatory, a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 
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as SN 1998bw are rare. Second, there is significant dispersion in the luminosities of Type 

Ib/c SNe, ranging from Lj^rad ~ 10^^ erg s~^ Hz^^ at the bright end (SN1998bw) to that of 
SN2002ap which is fainter by about four orders of magnitude (Berger, Kulkarni & Chevaher 
2002). It is curious that SN2002ap also happens to be the nearest Ib/c SN in our sample 
(Table 1). Six of the eight Type Ib/c SNe detected in the radio to date cluster in the range 
of about (3 — 50) x 10^^ erg s^^ Hz~^. This may be partly due to a selection effect since the 
typical detection threshold is about 4 x 10^^((i/50 Mpc)^ erg s~^ Hz"-*^. 

We also find that 28 of the 29 non-detections are no brighter than 0.1 times the lumi- 
nosity of SN1998bw. SN2002cg appears potentially brighter than SN 1998bw only because 
it is embedded in a radio bright host galaxy; we are therefore forced to use a lOcr limit on 
its luminosity (Table 1). Thus, the incidence of bright events like SN1998bw is < 3%. 

As with the radio luminosities, the peak times also exhibit great variation: at 1.4 GHz 
the emission from SN 2002ap peaked at about 7 days, while for SN 2002cj it peaked at about 
65 days. For SN 1998bw the initial peak occurred at 15 days, followed by a second peak at 
about 40 days. Similarly, at 8.5 GHz, SN 1998bw peaked at 12 and 30 days past explosion, 
SN1983N peaked at about 30 days, and SN2002ap is predicted to have peaked at ~ 1 — 2 
days (the first observation at this frequency was taken about 4 days after the SN explosion) . 



3.2. Expansion Velocities 

If the radio emission arises from a synchrotron spectrum peaking at the self-absorption 
frequency, f^, then the peak time and peak luminosity directly measure the mean expansion 
speed (Ghevalier 1998). This is simply because the self-absorption frequency is sensitive 
to the size of the source, while the luminosity is sensitive to the swept-up mass. We use 
Equation 16 of Ghevalier (1998) to evaluate the average expansion velocities: 

V, ^ 3.1 X lOX/f V.lo^pJ kms-^ (1) 

Here, Lp = 10^^Lp^26 erg s~^ Hz~^ is the peak luminosity, tp = lOt^^io days is the time of peak 
emission relative to the SN explosion, and Up — 5Up^^ GHz is the peak frequency. We infer 
velocities ranging from v ~ 10^ to 10^ km s~^ (Figure 3). Again, as with the luminosities, 
SN 1998bw with i> ~ c is an exception. 

We note that if free-free absorption plays a significant role, then i/p is only an upper limit 
to Ua, and Lp is a lower limit to the intrinsic peak luminosity. In this case, the inferred values 
of Vp listed above will in fact be a lower limit to the actual expansion velocity. However, 
this is probably not significant for Type Ib/c SNe since their compact progenitors have high 
escape velocities and therefore fast winds and low circumburst densities. Indeed, there is no 
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evidence for free-free absorption either for SN 1998bw (Kulkarni et al. 1998; Li & Chevalier 
1999) or SN 2002ap (Berger, Kulkarni & Chevalier 2002). 

Estimating the expansion velocities for the non-detections is not straightforward since we 
cannot ensure that the limits constrain the peak luminosity. We are therefore forced to make 
an additional assumption. For example, if we assume that most Type Ib/c SNe are similar 
in their emission properties to SN1983N, then the majority of upper limits approximately 
sample the peak emission and the inferred upper limits are < 0.3c (Figure 3). 

On the other hand, if the typical peak time is only a few days then our observations 
only constrain the decaying portion of the lightcurve and the expansion velocity may be 
higher. Fortunately, this is not a significant problem based on the following argument. The 
equipartition energy directly depends on the peak luminosity, f/cq ~ 3.7 x l^'^^L^^lj'^ erg 
(see §A.3), where Lp = 10^''Lp27 erg s~^ is the peak luminosity at 8.5 GHz. With a typical 
fading rate of Fj, oc in the optically thin regime, there are a few SNe that could have 
reached a peak luminosity of the order of 10^^ erg s~^ ii tp — 1 day post explosion. This 
is a reasonable minimum peak time taking into account the deceleration time of the ejecta. 
Thus, the equipartition energy is at most 10^^ erg, about two orders of magnitude lower than 
typical GRBs (§5). For most non-detections the hmit is in fact much lower, ~ 5 x 10^^ to 
2q48 gj,g ^]-^jg indicates that a few of the non-detected sources may have in fact produced 
mildly relativistic ejecta, but these would still be energetically uninteresting when compared 
to SN 1998bw let alone GRB afterglows. 



3.3. Energetics 

In the previous section we found that no SN observed to date is comparable to SN 1998bw 
especially in regard to the mean expansion speed. SN 1998bw was also interesting because 
it possessed an imusually large amount of energy in mildly relativistic ejecta. However, a 
purely hydrodynamic explosions can also produce some amount of relativistic ejecta. The 
energy of such ejecta can be estimated using well understood models of shock propagation 
in the pre-supernova cores (Chevalier 1982; Matzner & McKee 1999). The key parameters 
are Ek and Mej which can be inferred from the optical lightcurves and spectra using hydro- 
dynamic models of a SN explosion in a CO core coupled with radiative transfer calculations 
(e.g. Iwamoto et al. 2000). 

In this section, we investigate whether any of the detected Type Ib/c SNe possess such 
large energy in high velocity ejecta that cannot be explained by the simplest hypothesis 
of a purely hydrodynamic explosion. To this end, in Table 2 we summarize the results of 
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hydrodynamic models for the SNe that have been detected in the radio. 

The ejecta produced in a hydrodynamic explosion has a density profile that can be 
described by power laws at low and high velocities, separated by a break velocity, which for 
Type Ib/c progenitors is given by (Matzner & McKee 1999): 

^;ej,b ^ 5.1 X lQ'{EK,5i/M,j,iy/^-kms-\ (2) 

Here Ek = 10^^Ek,5i erg and M^j = 10Mej,i Mq. For typical values of Ek and Mej, the 
radio emission from the detected SNe is produced by ejecta above the break velocity. In 
particular, for SN2002ap, Vej,h ~ 2 x 10^ km s~^, which is lower than the velocity of the 
ejecta producing the radio emission, 1) 9 X 10^ km s"^ (Bcrgcr, Kulkarni & Chevaher 2002). 
Similarly, for SN 1998bw v^^h ranges from about 1.5 x 10"' to 3.5 x 10^ km (depending on 
which model is assumed. Table 2) while the radio emission was produced by ejecta expanding 
with r/3 PS 2. 

The ejecta velocity profile extends up to a cutoff determined by significant radiative 
losses when the shock front breaks out of the star. For a radiative stellar envelope this is 
^^cj,max ^ 11.5 X 10^^^^5\MtO-^2 (Matzucr & McKee 1999), assuming a stellar radius of 1 R©. 
For the SNe considered here we find cutoff velocities of F/? ~ 1 — 3. 

To determine whether there is sufficient energy in fast ejecta to account for the radio 
observations we calculate the energy above a velocity, V (Matzner & McKee 1999): 

E{v >V)^ -pfvHnvH^dv ^ 7.2 x lO^'^El^^M-Y'^V^'^-^^ erg, (3) 

where V5 is the velocity in units of 10^ km s~^. 

Unfortunately, as can be seen from Table 2, only four (including SN 1998bw) SNe have 
sufficient optical data which is necessary to estimate Ek and Mgj. Of this limited sample, 
much of the radio data for SN 19941 remain unpublished. Thus we are left with SN 2002ap, 
SN 1983N, and SN 1998bw. 

Using the parameters given in Table 2 for SN2002ap, Berger, Kulkarni & Chevalier 
(2002) find E{v > 0.3c) ~ 3.8 x 10"^^ erg. In contrast, from the radio observations we 
estimate 2 x 10^^ erg. Thus, there is no need, nor indeed room, for mildly relativistic ejecta 
in this SN. We therefore disagree with the claims of high velocity jets carrying a large amount 
of energy, ~ 0.1 - 1 FOE, made by Kawabata et al. (2002) and Totani (2003). Furthermore, 
the large discrepancy between the amount of energy inferred from the hydrodynamic models 
and the radio observations suggests that cither the optically-derived parameters arc in error. 
Equation 3 has an incorrect pre-f actor, or the radio estimate is incorrect. However, the radio 
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estimate is relatively robust (eventually related to equipartition energy estimates). On the 
other hand, as with SN 1998be (Hoflich, Wheeler & Wang 1999) the total kinetic energy may 
have been over-estimated, possibly as a result of neglecting a mild asymmetry. 

For SN 1984L we do not have a direct estimate of the energy in the radio-emitting ejecta 
since the peak of the radio emission has been missed. However, based on the similarity to 
SN1983N in the optically thin regime we estimate Lp{t = SOd.Up = 5 GHz) ^ 1.4 x 10^'^ 
erg s^^ Hz^^. This translates to a peak flux of 3.2 mJy at the distance of SN 1984L {d ^ 19 
Mpc). Using the equipartition analysis presented in §A.3 we estimate an energy of about 
7 X 10"^^ erg, and an average expansion velocity of about 0.1c. From Equation 3 we find 
E{v > 0.1c) ^ 3 X 10^° erg - similar to the conundrum discussed above for SN2002ap. 

For SN 1998bw, on the other hand, wc find E{v > c) ~ 2 x 10'^^ erg using the param- 
eters inferred by Hofiich, Wheeler & Wang (1999), or E{v > c) ^ 3 x 10*^^ erg using the 
parameters given by Iwamoto et al. (1998). In both cases, the energy available in fast ejecta 
is significantly lower than the energy inferred from the radio emission, ~ 10^° erg. 

To conclude, for SN 1984L and SN 2002ap a hydrodynamic explosion can supply the 
energy and velocity that are responsible for the observed radio emission. Most likely, the 
same is true for the non-detections. On the other hand, SN 1998bw is a clear exception, 
exhibiting a significant excess of energy in ejecta moving with F/? 2 compared to what is 
available from a hydrodynamic explosion. 

4. A Compcirison to 7-Ray Burst Afterglows 

In the previous section wc investigated the radio properties of Type Ib/c SNc and found 
that in every respect SN 1998bw was unique. In this section we compare the Ib/c sample 
(including SN 1998bw) with the radio afterglows of GRBs. In Figure 2 we plot the radio 
hghtcurves of GRB 970508 (Frail, Waxman & Kulkarni 2000), a typical cosmological burst, 
and the nearest event, GRB 030329 (Berger et al. 2003), in addition to the SN hghtcurves. 
As demonstrated by this figure and Figure 4, the radio hghtcurves of GRB afterglows and 
SNe are dramatically different. Furthermore, SN 1998bw is unique in both samples: it is 
fainter than typical radio afterglows of GRBs but much brighter than Type Ib/c SNe. 

Figures 2 and 4 have significant imphcations, namely none of the Type Ib/c SNe pre- 
sented here could have given rise to a typical 7-ray burst. It has been suggested that GRBs 
are distant Type Ib/c SNe but with their jets pointed at the observer, whereas such a bias 
is absent in the nearby Type Ib/c sample. However, most of our radio observations are 
obtained on a timescale of 10-100 days (see Figure 1). Scaling from the observed "jet" break 



-10- 



times of a few days in GRB afterglows, off-axis coUimated explosions become spherical on 
a timescale of ~ 10 — 10^ days (Paczynski 2001; Granot & Loeb 2003) at which point the 
relative geometry between the observer and the explosion is not important. Thus, we find 
no evidence suggesting that all or even a reasonable majority of Type Ib/c SNe give rise to 
GRBs. We now quantify the difference between the Ib/c and GRB samples. 

Our goal here is to determine the differential luminosity distribution, n{L), which agrees 
with both detections and upper limits; here n{L) is the number of events with luminosity 
between L and L + dL. It is important to include non-detections since they represent the 
majority of the SN data. Similarly, we include upper limits on the radio luminosity of GRB 
afterglows that have been localized in other wave-bands (i.e. optical and X-rays) and for 
which a redshift has been measured. Unfortunately, as many as half of the GRBs localized 
in the X-rays do not have a precise position, and hence a redshift. For these afterglows it is 
not possible to provide a limit on the radio luminosity. Still, with a typical flux limit of about 
0.3 mJy (5cr; Frail et al. 2003), and assuming that these sources have a similar distribution 
of redshifts to the detected afterglows, we find typical luminosity limits of about 10^^ erg 
s~^ Hz^^, consistent with the peak of the distribution of detected afterglows. Therefore, 
unless these sources are biased to low redshift we do not expect a strong bias as a result of 
neglecting them. 



The quality of fit for n(L) is determined using the Likelihood function, C — ni=i 
with (Reichart & Yost 2001): 

. _ / J^^n{L)Q{Li,aLi)dL = detection 



where N is the total number of sources (SNe or afterglows), and G{Li,aLi) is a normalized 
Gaussian profile centered on the observed luminosity of a detected source and with a width 
equal to the la rms uncertainty in the luminosity. 

We consider four models for n{L) based on the apparent distribution of the detections 
and upper limits: a Gaussian, 
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a decreasing power-law. 




L< Lo 
L>Lo 
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an increasing power-law. 




L"V(^o^ ) < L < L, 
L>Lo 



(7) 



and a flat distribution, 



r L < Li 

n{L) = I 1/(L2 - Li) Li<L<L2 . (8) 
i L > L2 

In each case we fit for tlic two free parameters (e.g. Lq and ctl in Equation 5). We do not use 
the increasing power law model for the individual distributions since the observations are 
clearly inconsistent with such a model. The resulting best-fit models are shown in Figure 5 
and summarized in Table 3. 

We find that the SN population is modeled equally well with the Gaussian, fiat, or 
decreasing power law distributions, while the GRB afterglows can be fit with a Gaussian 
or fiat distributions; a decreasing power law provides a much poorer fit. Regardless of the 
exact distribution the two populations require distinctly different parameters, with a minimal 
overlap at the tails of the distributions. 

Fitting the combined SN and afterglow data with the models provided above (Figure 5 
and Table 3) we find that even the best models (an increasing power-law or a fiat distribu- 
tion) provide a much poorer fit; the likelihood of the fits are ln(>C) ~ 104 compared to the 
combined value of 61.3 for the separate Gaussian fits. Thus, the two populations cannot be 
accommodated with a simple single distribution. This points to a separate origin for the 
GRB and Type Ib/c SN populations. However, SN 1998bw can be accommodated in either 
population. It is equally plausibe that it is a low luminosity GRB or the brightest Type Ib/c 
radio supernova known to date. 



5. Discussion and Conclusions 

We presented VLA radio observations of 33 Type Ib/c SNe observed between late 1999 
and the end of 2002. Four of these SNe have been detected, giving a detection rate of about 
12% above a typical 3a fiux limit of 0.15 mJy. At the same time, the combined detections 
and non-detections indicate that at most 3% of Type Ib/c SNe are as luminous as SN 1998bw, 
although the single source which may be brighter is only so because it is embedded in a radio 
bright host galaxy. 

We infer typical velocities of the radio-emitting ejecta of about 10^ — 10^ km s~^ for 
the detected SNe, and upper hmits in the same range for the non-detections. We also find 
that a hydrodynamic explosion can supply the energy carried by the fastest ejecta. Finally, 
none of the detected SNe show clear evidence for variable energy input (shells with different 
velocity or continued activity by the central engine); however, we note that our sampling is 
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quite sparse. 

The measurements (radio lightcurves) and inferences (energy in fast ejecta, energy addi- 
tion) offer no compelling reason to conclude that any of our SNe have the special properties 
of SN1998bw (§5.1). 

Norris (2002) has proposed, based on the empirical lag-luminosity relation, that 25% of 
Type Ib/c SNe are similar to SN 1998bw. 

We also compared the Type Ib/c sample with the sample of radio afterglow of GRBs. 
Empirically, these two populations appear to be quite disparate. This conclusion is reinforced 
by careful modeling of the luminosity distributions. Still, SN 1998bw may belong to either 
population. 

5.1. What is SN 1998bw? 

Our four year survey of Type Ib/c SNc was first and foremost motivated by the peculiar 
object, SN 1998bw. This supernova showed three attributes unique to GRBs: relativistic 
ejecta, substantial reservoir of energy in such ejecta, and energy addition. A singular hy- 
drodynamic explosion cannot account for these attributes. A natural explanation is that 
SN 1998bw, like GRBs, was , driven by an engine powerful enough to significantly modify the 
explosion. 

Our survey has demonstrated that SN 1998bw-like events are rare in the local sample. 
This begs the question: what is SN 1998bw? 

Two popular scenarios have been suggested. The first - the "off-axis" scenario - holds 
that SN 1998bw is a typical GRB albeit nearby and with coUimatcd ejecta pointed away from 
us (Hoflich, Wheeler & Wang 1999; MacFadyen & Woosley 1999; Nakamura 1999; Granot 
et al. 2002). This hypothesis is attractive because of its simphcity. We know GRBs exist 
and most of them do not point towards us (Prail et al. 2001). 

In the other scenario SN 1998bw is a new type of explosion with little energy in ultra- 
relativistic ejecta (Bloom et al. 1998; Kulkarni et al. 1998). Evidence in favor of this idea 
is best illustrated by Figure 6 where we find that GRB 980425 is consistently at the faint 
end of the GRB population. 

Unfortunately, we are not able to decisively resolve this controversy. As demonstrated 
by Figure 4, one could argue that SN 1998bw is at the bright end of the radio luminosity 
function of Type Ib/c supernovae or at the faint end of GRB radio afterglow. 
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The expected rate of SN 1998bw-like events in the off-axis framework is about 0.5% of 

Type Ib/c events, given the average beaming factor of about 500 derived by Frail et al. 
(2001). Thus, it is not entirely improbable that one out of about 40 Type Ib/c SNe observed 
to date is an off-axis GRB. As an aside, we can use our 3% limit and compare the event rate 
of Type Ib/c SNe with the observed rate of GRBs (§1) to place a limit of /{, > 3 x 10"'^ on 
the beaming fraction. This corresponds to a limit of 9j > 1.4° on the jet opening angles of 
GRBs; narrower jets are not likely. This result may also be interpreted as a limit on angular 
size of the highly relativistic core in models of variable energy and/or Lorentz factor across 
the surface of the jet (Rossi, Lazzati & Rees 2002; Zhang & Meszaros 2002). 

In the gSN framework, we now know that at most a few percent of Type Ib/c are 
possibly gSNe. At the same time, the recent GRB 030329 was accompanied by a SN similar 
to SN1998bw (SN2003dh; Stanek et al. 2003; Hjorth & et al. 2003). Thus, an investigation 
of the number and properties of gSNe requires observations of both local Type Ib/c SNe and 
GRBs. 

While we cannot determine the exact origin of SN 1998bw based on the statistics of our 
survey, the ultimate detection of similar events at the level of about 1% may in fact allow 
us to distinguish between the off-axis and gSN scenarios. 

5.2. Hypernovae 

The discovery of broad optical lines in SN 1998bw and large explosive energy release, 
> few FOE, prompted some astronomers to use the designation "hypernovae" for SN 1998bw- 
hke SNe. Unfortunately, this designation is not well defined. To begin with, the term 
hypernova was first used by Paczynski (1998) to describe the GRB/afterglow phenomenon; 
thus, this term implies a connection to GRBs. The prevalent view now is that hypernovae 
are characterized by broad optical absorption lines and larger than normal energy release. 
However, neither of these criteria has been defined quantitatively by their proponents. 

Ignoring this important issue, the following have been suggested to be hypernovae: the 

Type Ib/c SNe 1992ar (Glocchiatti et al. 2000), 1997dq (Mathcson et al. 2001), 1997ef 
(Iwamoto et al. 2000; Mazzali, Iwamoto & Nomoto 2000), 1998ey (Garnavich et al. 1998), 
and 2002ap (Mazzah et al. 2002), and the Type II SNe 1992am (Hamuy 2003), 1997cy 
(e.g. Germany et al. 2000), and 1999E (Rigon et al. 2003). Some have also been claimed 
to be associated with GRBs detected by BATSE, but at a low significance. 

Our view is that the critical distinction between an ordinary supernova and a GRB 
explosion are relativistic ejecta carrying a considerable amount of energy. Such ejecta are 
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simply not traced by optical spectroscopy. This reasoning is best supported by the fact 
that the energy carried by the fast ejecta in SN1998bw and SN2002ap (Bergcr, Kulkarni 
& Chevaher 2002) differ by four orders of magnitude even though both exhibit broad spec- 
tral features at early times. Thus, broad lines do not appear to be a good surrogate for 
SN 1998bw-like objects. 

In addition, in two cases, SNe 2002ap and 1984L, the energy inferred from the radio 
observations indicates that the total kinetic energy as inferred by optical spectroscopy and 
light curves may have been over-estimated by an order of magnitude (§3.3). It is possible 
that "hypcrnovac" arc in fact only slightly more energetic than typical Type Ib/c SNe, but 
exhibit a mild degree of asymmetry, leading to excessively high estimates of the total energy 
when a spherical explosion is assumed. 

We suggest that the term hypernova be reserved for those SNe, like SN 1998bw, which 
show direct evidence for an engine through the presence of relativistic ejecta. As illustrated 
by SN 1998bw, the relativistic ejecta are reliably traced by radio observations. 

We end with the following conclusions. First, radio observations provide a robust way 
of measuring the quantity of energy associated with high velocity ejecta. This allows us 
to clearly discriminate between engine-driven SNe such as SN 1998bw and ordinary SNe, 
powered by a hydrodynamic explosion, such as SN2002ap (Berger, Kulkarni & Chevalier 
2002) and SNe 2001B, 2001ci, and 2002cj presented here. Second, since at least 97% of local 
Type Ib/c SNe are not powered by inner engines and furthermore have a total explosive 
yield of only 10^^ erg in fast ejecta (where fast means v ~ 0.3c compared to F ~ few 
in GRBs), there is a clear dichotomy between Type Ib/c SNe and cosmic, engine-driven 
explosions (Figure 6). The existence of intermediate classes of explosions and the nature of 
SN 1998bw may be ascertained with continued monitoring of several hundred Type Ib/c SNe 
and cosmological explosions. Fortunately, such samples will likely become available over the 
next few years. 

We thank B. Schmidt for helpful comments. GRB and SN research at Caltech is sup- 
ported in part by funds from NSF and NASA. 
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A. Results for Individual Supernovae 
A.l. SN2001B 

SN 2001B was discovered in images taken on 2001, Jan 3.61 and 4.57 UT, approximately 
5.6 arcsec west and 8.9 arcsec south of the nucleus of IC391 (Xu & Qiu 2001). The SN 
explosion occurred between 2000, Dec 24.54 UT and the epoch of discovery. Based on an 
initial spectrum, taken on 2001, Jan 14.18 UT Matheson et al. (2001) concluded that the SN 
was of Type la approximately 9 days past maximum brightness, showing well defined Si 11 
and Call features, with a Sill expansion velocity of about 7400 km sec^^. A subsequent 
spectrum taken on 2001, Jan 23 indicated that the SN was in fact of Type lb based on clear 
Hel absorption lines (Chornock & Filippenko 2001). 

We observed SN2001B on 2001, Feb 4.49 UT at 8.46 GHz. An initial analysis revealed a 
source which was interpreted as the host galaxy of the SN, with a flux of 3.5 mJy A second 
epoch obtained on 2002, Oct 28.45 UT revealed that the source has faded below 0.12 mJy, 
estabhshing it as the radio counterpart of SN2001B. (Figure 1). 

A.2. SN2001ci 

SN2001ci was initially detected in images taken with the Katzman Automatic Imaging 
Telescope on 2001 Apr 25.2 UT, 6.3 arcsec west and 25.4 arcsec north of NGC3079 (Swift, Li 
& Filippenko 2001). These observations did not provide conclusive evidence that the source 
was in fact a SN. A spectrum obtained by Filippenko & Chornock (2001) on 2001, May 30 
UT revealed that the source was in fact a Type Ic SN about 2 — 3 weeks past maximum 
brightness (Matheson et al. 2001). 

We observed the SN on 2001, Jun 6.35 UT at 8.46 GHz, but did not detect the source 
since it appeared to be part of the host galaxy extended structure. A second epoch taken 
on 2002, Jun 10.9 UT revealed a clear fading at the optical position of the SN with a flux of 
1.45 ± 0.25 mJy in the flrst epoch, and a 2a hmit of 0.3 mJy in the second epoch. 

A.3. SN2002cj 

SN2002cj was discovered on 2002, Apr 21.5 UT, 1.4 arcsec west and 3.9 arcsec south 
of the nucleus of ESQ 582-G5, at a distance of 106 Mpc (Ganeshahngam & Li 2002). The 
SN explosion occurred between Apr 9.5 UT and the epoch of discovery. Spectra of the SN 
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obtained on 2002, May 2.43 and May 7 UT revealed that SN2002cj is of lype Ic (Matheson 
et al. 2002; Chornock 2002). 

We initially observed the SN on 2002, Jun 3.19 UT at 1.43, 4.86, and 8.46 GHz and 
detected a faint source at all three frequencies. The SN spectrum is given by F^, oc i/-o i±o-3 
between 1.43 and 4.86 GHz, and oc i/-i-2±o.8 between 4.86 and 8.46 GHz, indicating that 

the spectral peak was most likely located between 1.43 and 4.86 GHz during the first epoch 
(At 43 — 55 days). In subsequent observations at 1.43 GHz the SN brightened and then 
faded, as expected if the peak was in fact above 1.43 GHz initially, and shifted through the 
band at later epochs. Using the expected shape of the spectrum, F^, oc z/^/^ for u < Up 
and Fj, oc for u > Up (with p ~ 3), we find F^^p ~ 0.5 mJy at z/p ~ 2 GHz and 

= 43 - 55 days. 

We use these values along with the well-established equipartition analysis (Readhead 
1994) to derive some general constraints on the properties of the emitting material. In 
particular, the energy of a synchrotron source with flux density, Fp(i/p, tp), can be expressed 
in terms of the equipartition energy density, 

^ = lesv'\l + ^V-''), (Al) 

Wherer^ = ^,/^eq, ^eq «^ 120dMpf i^£j>S"''^^'' /^a^' = l-lXl05«d^p,F^^,^jy^-^HzC.a. 

erg, and eg and sb are the fractions of energy in the electrons and magnetic fields, respectively. 
In equipartition eg = = 1 and the energy is minimized; a deviation from equipartition 
will increase the energy significantly. 

For SN2002cj we find e^q{t = 43 - 55 d) fa 30 //as (i.e. R^^ 5 x 10^^ cm), which 
indicates an average expansion velocity, fcq ~ (0.35 — 0.45)c. The equipartition energy is 
Ueq 8 X 10^^ erg, indicating a magnetic field strength of Beq ~ 0.2 G. 
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Table 1. Radio Observations of Type Ib/c Supernovae in the Period 1999-2002 
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Nov 18.02 
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54 


No 
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±71 
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59 
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2000S 
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<Feb 28^ 
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7443 


Jun 21.2-25.9 
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15 
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72 
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3500 ± 29 
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269.9-469 
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±85 


2001em 


7722 


Sep 10.3-15.3 
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Table 1 — Continued 
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Dist. 


Detected? 
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2001fx 
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7800 


Jan 9.4-21.4 


Jun 18.0 
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58 
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Feb 1.03 


3.5 


7 


Yes 


374 ± 29 
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±66 
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Jun 3.19 


42.7-54.7 
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112 ±33 
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57.5-69.5 
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Jun 26.0 


65.5-77.5 




Yes 






300 ± 68 


2002cp 
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No 


±31 
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2002CW 
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2001 Oct 1.2 May 16.5 






71 










2002dg 
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2002 May 5-May 31.3 


Jul 9.95 


39.7 66 


215' 
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±33 


92 ±37 


±46 


2002dn(?) 
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May 31.5-Jun 15.5 


Jul 4.4 


18.9-33.9 
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No 


±39 


±43 


±72 


2002dz 


7935 


2001 Nov 10.2-Jul 16.5 






84 










2002ex 
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Aug 19.3 
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2002fli(?) 
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Apr 9-May 9.3 
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2002ge(?)'" 
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~ 13 


47 


No 
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2002gy 
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Oct 9.4-16.4 


Oct 28.45 


19.1-26.1 
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±36 
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2002hf 
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Oct 22.3-29.3 


Nov 7.17 
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No 


±40 






2002hn 
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Oct 21.5-30.5 


Nov 7.42 
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No 


±44 






2002ho 
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May 27 Nov 5.1 


Nov 15.33 


~ 56P 


42 


No 
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2002hy 
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Nov 21.71 


9.6-39.6 


58 


No 


±74 






2002hz 
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Nov 2.2-12.2 


2003 Jan 21.09 


69.9-70.9 


85 


No 


±29 






2002ji 
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Apr 10-Nov 30.8 


Dec 5.67 


38.4-59.4« 


23 


No 


±43 






2002jj 
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Oct 1.4-24.4 


Dec 15.28 


51.9-74.9 


66 


No 


±33 






2002jp 
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May 14.2-Nov 23.5 


Dec 14.55 


~ 33'^ 


58 


No 
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SN 


lAUC 


to 


^obs 


At 


Dist. 


Detected? 


-^8.46 


^4.86 


^1.43 






(UT) 


(UT) 


(days) 


(Mpc) 








(/^Jy) 


2002jz 


8037 


2001 Dec 5-2002 Dec 23.3 


2003 Jan 3.28 


~ 32^ 


24 


No 


±25* 







Note. — The columns are (left to right): (1) SN name, (2) lAU Circular number for the initial detection, (3) time 
of the SN explosion, with the range given between the most recent observation of the galaxy which did not show the 

SN and the epoch at which the SN was detected, (4) epoch of our VLA observations, (5) time delay between the SN 
explosion and the epoch of our observations, (6) distance to the galaxy (assuming Hq = 65 km s~^ Mpc~^), (7) indicates 
whether radio emission was detected, (8) flux density at 8.46 GHz, (9) flux density at 4.86 GHz, and (10) flux density 
at 1.43 GHz. 

° uncertainties are quoted as la rms; nebular phase; ^ falls on top of galaxy substructure, < lOcr; la or Ic (lAUC 
7574); ^ before maximum (lAUC 7574); ^ falls on top of galaxy with complex substructure, < 4(7; ^ falls on top of 
galaxy, < 5a; ^ falls on top of galaxy, < 35f7; ' on top of galax;y substructure, < 3a; ^ See Berger, Kulkarni & Chevalier 
(2002); on top of galaxy, < lOcr; ' lAUC 7922; ™ Type Ic similar to SN 19941 a few days before maximum, or a 
sub-luminous Type la (lAUC 7990); " on top of galaxy, < 10a; ° on top of galaxy, < 3a; ^ within a few weeks past 
maximum brightness (lAUC 8014); '3 — 6 weeks past maximum light (lAUC 8026); ^ two weeks past maximum light 
(lAUC 8031); * ten days past maximum (lAUC 8037); * on top of galaxy substructure, < 3a. 
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Table 2. Ejecta and Progenitor Properties of Type Ib/c Supernovae Detected in the Radio 
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jTvei. 




l-LU ergj 












1984L 


20 


50 


0.2 






1 






10 








2 








0.1 




20-30 


3 


19941 


1 - 1.4 


0.9-1.3 


0.07 


< 1.5 


< 15 


4 






0.9 


0.07 


1.35 




5 




1 










6 


1998bw 


30 




0.7 


13.8 


40 


7 




50 


10 


0.4 


13.8 


40 


8 




22 




0.5 


6.5 




9 




2 


2 


0.2 






10 








0.5-0.9" 






11 


2002ap 


4-10 


2.5-5 


0.07 


5 


20-25 


12 



Note. — The columns are (left to right): (1) Ejecta kinetic energy, 
(2) ejecta mass, (3) ^^Ni mass, (4) estimated mass of the CO core, 
(5) estimated mass of the progenitor, and (6) references. Data are 
not available for the SNe detected in this survey. " These authors 
use the models of Iwamoto et al. (1998) and Woosley, Eastman & 
Schmidt (1999) as input; they assert that observations in the nebular 
phase exclude the low ^^Ni mass inferred by Hofiich, Wheeler & Wang 
(1999). 

References. — (1) Baron, Young & Branch (1993); (2) Swartz & 
Wheeler (1991); (3) Schlegel & Kirshner (1989); (4) Young, Baron & 
Branch (1995); (5) Iwamoto et al. (1994); (6) Millard et al. (1999); 
(7) Iwamoto et al. (1998); (8) Nakamura et al. (2001); (9) Woosley, 
Eastman & Schmidt (1999); (10) Hofiich, Wheeler & Wang (1999); 
(11) SoUerman et al. (2000); (12) Mazzah et al. (2002) 
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Table 3. Best-Fit Models for the Supernova and 7-Ray Burst Populations 



Population 


Model 


Parameters 


log(>C)/dof 


SN 


Gaussian 


(26.1, 1.0) 


22.5/36 


SN 


D. Power law 


(-29.0, 25.4) 


22.4/36 


SN 


Flat 


(20.0, 29.1) 


22.3/36 


GRB 


Gaussian 


(31.0, 0.8) 


38.8/33 


GRB 


D. Power law 


(-22.5, 29.6) 


48.9/33 


GRB 


Flat 


(29.6, 32.4) 


37.7/33 


SN+GRB 


Gaussian 


(28.4, 2.4) 


130.3/71 


SN+GRB 


D. Powerlaw 


(-10.4, 25.3) 


123.8/71 


SN+GRB 


1. Powerlaw 


(3.7, 32.3) 


103.7/71 


SN+GRB 


Flat 


(22.6, 32.4) 


104.9/71 



Note. — The columns are (left to right): (1) Data set, 
(2) population distribution function, (3) best-fit parame- 
ters, and (4) log hkelihood. A detailed explanation of the 
models and the fitting procedure is provided in §4. 
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Fig. 1. — Radio lightcurves of Type Ib/c SNc detected in this survey and from the hterature, 
as well as upper limits for the non-detections; these are plotted as 3a in most cases with 
the exception of SNe which are located on top of a bright host galaxy (see Table 1). For 
SN2002ap we plot the model of Berger, Kulkarni & Chevaher (2002), while the other solid 
lines simply trace the observations and do not represent a model fit. The uncertainty in 
time for the non-detections represents the uncertain time of explosion. We note that for 
SN2002cg, which is the only SN that is potentially brighter than SN1998bw, the hmit is 
10(7 due to the superposition of the SN on top of its host galaxy. 
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Fig. 2. — Same as Figure 1 but including the 8.46 GHz lightcurves of GRB 970508 (Frail, 
Waxman & Kulkarni 2000) and GRB 030329 (Berger et al. 2003). These GRB afterglows 
are at least two orders of magnitude brighter than SN 1998bw, the brightest Type Ib/c SN. 
We note that the fluctuations in the GRB lightcurves are not intrinsic, and arise instead 
from interstellar scintillation. Based on the significant difference in radio luminosity we rule 
out the possibility that the Type Ib/c SN observed here produced a GRB. This is discussed 
more quantitatively in in §4 and Figures 4 and 5. 
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Fig. 3. — Peak radio luminosity plotted versus the time of peak luminosity for Ib/c SNe 
studied in this survey and from the literature. Symbols are as in Figure 1. The diagonal 
lines are contours of constant average expansion velocity based on the assumption that the 
peak of the radio luminosity occurs at the synchrotron self- absorption frequency (Chevalier 
1998). While the upper limits do not necessarily measure the peak of the spectrum at the 
time of the observation, a comparison to SN 1983N indicates that the range of time delays 
relative to the SN explosion reasonably samples the peak. Upper limits measured at t > 100 
days probably miss the peak of the synchrotron spectrum and therefore do not provide a 
useful limit. 
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Fig. 4. — Histograms of the radio luminosity of Type Ib/c SNe from this survey and the 
hterature, and GRB radio afterglows from the sample of Frail et al. (2003). Upper limits are 
plotted as 3a, unless there is significant contamination from the host galaxy (see Table 1). 
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Fig. 5. — Same as Figure 4 but with models of the luminosity distribution (§4). The dashed 
lines are a Gaussian profile fit for the detections in each sample. The solid lines are a fit to 
the detections and upper limits using several models for the distribution function (§4). In 
the top panel we model each population separately, whereas the bottom panel shows models 
for the combined populations. No single distribution can fit both the local Type Ib/c SN 
population and the cosmological GRB population (Table 3). 
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Fig. 6. — Histograms of the beaming-corrected 7-ray energy (Bloom, Frail & Kulkarni 2003), 
E'-y, the kinetic energy inferred from X-rays at i = 10 hr (Berger, Kulkarni & Frail 2003), 
Ek,x, and total relativistic energy, + Ek-, where Ek is the beaming-corrected kinetic 
energy inferred from the broad-band afterglows of GRBs (Li & Chevalier 1999; Panaitescu & 
Kumar 2002) and radio observations of SNe. The wider dispersion in E"^, and Ek,x compared 
to the total energy indicates that engines in cosmic explosions produce approximately the 
same quantity of energy, thus pointing to a common origin, but the ultra-relativistic output 
of these engines varies widely. In Type Ib/c SNe, on the other hand, the total explosive yield 
in fast ejecta (typically ~ 0.3c) is significantly lower. This points to a separate origin for 
these two explosive phenomena. 



